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The radio-frequency quadrupole linear ion trap is a widely used device in physics and chemistry. When used for trapping 
of large ion clouds, the presence of anharmonic terms in the radiofrequency potential limits the total number of stored 
ions. In this paper, we have studied the anharmonic content of the trapping potential for different implementations of a 
quadrupole trap, searching for the geometry best suited for the trapping of large ion clouds. This is done by calculating 
the potential of a real trap using SIMION8.0, followed by a fit, which allows us to obtain the evolution of anharmonic 
terms for a large part of the inner volume of the trap. 
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1. Introduction 

The trapping of charged particles by radiofrequency 
(RF) electric fields was first demonstrated in 1954 [1[ and 
quickly proved to be an extremely powerful tool for the 
experimental investigation of a wide range of phenomena. 
In particular the linear quadrupole RF trap is found at 
the heart of many experiments where few to many ions are 
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laser-cooled to very low temperature for applications in O i 
tical frequency metrology 0] , quantum computation 0. 
and formation of large Coulomb crystals [5]. In this paper, 
the deviation from the ideal harmonic potential is evalu- 
ated for different implementations of linear quadrupole ion 
traps, in order to maximize the number of trapped ions. 

In this article, we briefly recall the basic information 
about linear quadrupole traps; detailed descriptions of RF 
traps and trajectories of the stored ions can be found in 
several textbooks (e.g. The most simple linear RF 

trap is formed by four parallel rods of radius r e , see fig- 
ure [TJ In such devices, the trapping of charged particles in 
the transverse direction is achieved by applying an alter- 
nating potential difference between two pairs of electrodes 
(called the RF electrodes in the following). The trapping 
along the symmetry axis of the trap (to be taken as the z 
axis) is realised by a static potential applied to two elec- 
trodes sets at both ends of the trap (called the DC elec- 
trodes in the following). The potential at the centre of a 
quadrupole linear trap is then well approximated by: 




Figure 1: Schematic transverse view of a linear RF trap 
made of four cylindrical electrodes. The axes convention is 
shown, as well as the definitions of the applied RF voltages 
V r f = c/o+v ° 2 cosnt , the trap radius r and the electrode 
radius r e . 
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where Vq is the amplitude of the RF potential difference 
between two neighbouring electrodes, Uq a static contribu- 
tion to this potential and Udc the voltage applied to the DC 
electrodes. SI is the angular frequency of the RF voltage, ro 
is the inner radius of the trap, 2zo is the length of the trap, 
corresponding to the distance between the DC electrodes, 
k and C are geometric loss factors, k is mainly induced 
by shielding effects whereas C stands for a loss in trapping 
efficiency due to the non-ideality of the quadrupole Q- 
While in most cases <ft is an accurate approximation close 
to the trap axis and around z ~ 0, strong deviations oc- 
cur further away from the center. The deviations induced 
depend strongly on the geometry used for the trap, as it 
is shown in the present article. 

The motion of a single ion in a potential defined by 
4>(x, yA) is described by the Mathieu differential equa- 
tions |8J. The stability of the trajectory is controlled by 
the Mathieu parameters, which depend on the working pa- 
rameters of the trap and on the charge (Q) and mass (M) 
of the particle: 



q x = -q y 

a x = a — Aa 
4QU 



2QV 
MSl 2 r 2 C 



a 



MSl 2 r 2 C 



Aa 



-(a + Aa) 
~~ MzlSl 2 



(2) 



where a represents the Mathieu parameter of the pure 
quadrupole mass filter, and Aa is the perturbation intro- 
duced by the presence of the axial confining voltage 
Most traps are operated in the adiabatic regime, where 
it is possible to decouple the short time scale of the RF 
driven motion, called micromotion, and the longer time 
scale of the motion induced by the envelope of the elec- 
tric field, known as macromotion [loj . In this frame, the 
ion's motion is the superposition of a periodic oscillation 
in an effective harmonic potential, characterized by the 
radial frequencies w x (a x , q x , SI) and uj y (a yj q y ,Sl), with the 
RF-driven motion (at frequency SI) with an amplitude pro- 
portional to the local RF electric field. 

The simultaneous trapping of two or more ions cou- 
ples the x and y equations of motion through a nonlinear 
term induced by the Coulomb interaction and an analyt- 
ical solution does not exist any longer. Moreover, due 
to this nonlinear contribution the ions gain energy from 
the RF field. This energy gain depends on the ion den- 
sity and / or neutral background pressure inducing colli- 
sions While it can be negligible for very dilute ion 
clouds in good vacuum conditions, it can be a critical point 
for the storage of dense clouds where an important fraction 
of the ions can be lost due to this gain of energy. 



Nevertheless, even in the case of a single ion in a perfect 
vacuum, nonlinearities can be introduced in the equations 
of motion by differences between the real RF and DC po- 
tentials and the ideal case. Imperfections in the potentials 
come from misalignments of the different electrodes and 
differences between the electrode surfaces from the ideal 
cquipotential surfaces. These effects are very well known 
for hyperbolic Paul traps where holes in the electrodes are 
necessary for the injection and extraction of ions, electrons 
or light. Moreover, these undesired contributions to the 
potential can contain higher harmonic contributions which 
lead to a special form of RF-heating. More precisely, it has 
been shown I2J that for an imperfect Paul trap, nonlinear 
resonances leading to the loss of ions occur in the other- 
wise stable region, if the ion's secular frequencies and the 
radiofrequency are related by a linear combination with 
integer coefficients (N x ,N z ,k): 



N x uj x + N z uj z 



kSl. 
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This resonant phenomena due to nonlinear couplin gs h as 
been observed experimentally by several groups [UH]. 

Ideal harmonic potentials are obtained only if the shape 
of the RF electrodes follows the theoretical equipotential 
surface, a hyperbolic surface in the case of a quadrupole 
trap. The difficulty in machining such surfaces can be over- 
come by realizing circular rods. In this case, the relation 
r e = f.I45IIro between the trap radius ro and the elec- 
trode radius r e guarantees a minimal contribution of the 
lowest order term in the perturbation potential expansion 
for infinitely long RF electrodes This radius ratio is 
used for all the traps studied in this article. While some 
work has been done in optimizing miniature Q and cylin- 
drical [lU traps designs, to our knowledge no study has 
concentrated on the study of the anharmonic terms in a 
real linear trap. 

Our research project aims at trapping a very large and 
laser-cooled ion cloud to study fundamental phenomena 
related to their dynamics and thermodynamics. In order 
to reach a very high temporal stability of the ion number, 
the trapping potential has to be quasi-ideal over a large 
part of the trapping volume. The ion cloud is expected 
to fill the trap to half the trap size in both directions 
(— zq/2 < z < zo/2 and < r < rp/2) with a sufficiently 
high ion density such that a high number of ions (~ 10 7 ) 
can be reached in a compact trap. 

The present study evaluates the anharmonic contribu- 
tions to the trapping potential of various existing devices 
and proposes an optimised alternative version. Contribu- 
tions due to the finite size of the RF electrodes and to the 
shape of the DC electrodes are quantified and compared. 
Contrary to former studies in spherical traps 0] where an- 
harmonic contributions were evaluated in the very centre 
of the trap (r/ro,z/zo < 0.1), which is relevant for single 
ions or chains of single ions, the present study is made in 
three dimensions and extends to values beyond zq/2. This 
is motivated because a large ion cloud explores a much 
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Figure 2: (Color Online) Design of the first two analysed 
geometries. The dimensions used in both cases are: vq = 4 
mm and zq = 10 mm, with r e = 1.14511ro. Trap #1 is 
formed by four segmented rods whose electrical connec- 
tions are detailed in table [TJ Trap #2 consists of four con- 
tinuous rods surrounded by DC electrodes made of two 
external rings of longitudinal dimension l T i ng = 4 mm, 
inner radius rj„t = (tq + 2r e ) + 1 mm and outer radius 
Text = r int + 5 mm. The distance between these rings de- 
fines 2zo- The dotted lines define the volume used for the 
SIMION8.0 calculations. 

Table 1: Electrical configuration for the different trap ge- 
ometries. 
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bigger volume of the trap, and therefore it is desirable 
that this volume has a low anharmonic content in order to 
minimise RF-heating. 

Figure [2] shows the two geometries initially analysed. 
We have chosen only devices which leave the z-axis free as 
this axis is often needed to implement laser cooling and/or 
to introduce/extract ions. The segmented rod geometry 
(#1) is used, among others, by the Ion Trap Group in 
Aarhus University to trap large ion crystals [a, |9( . 

This geometry can be used with two different electri- 
cal configurations, named #1.1 and #1.2 (see table [1] for 
details). #1.2 has been analyzed due to its simplicity as 
#1.1 requires a much more advanced electronic configura- 
tion than #1.2 (see tabled]). The second implementation, 
#2, is similar to one of the trap used in Quantum Op- 
tics and Spectroscopy Group at Innsbruck University for 
quantum information experiments 17 j. 

This article is organised in several sections. The method 
used to calculate the potential created by a given geom- 
etry and the fitting procedure is described in section [2] 
In section [3l the anharmonic contributions are calculated 



and evaluated for the RF-part while we discuss the DC 
potential in section Finally, an alternative geometry 
is introduced and analysed in section [5j followed by the 
conclusion, section [S] 

2. Fitting the "real" trap potentials 

2.1. Computation of the potentials 

The commercial software SIMION8.0 18] has been used 
to numerically solve the Laplace equation for each geome- 
try. SIMION8.0 uses a Finite Difference Method (FDM), 
where the solution is obtained at each node of the mesh 
used to describe the electrodes and the space between 
them. The error on the electrostatic potential using FDM 
scales with l/h [l^|, where h is the node density (number 
of nodes/mm), but the reached accuracy depends on the 
particular geometry of the electrodes. In order to check the 
l/h dependency and to obtain an estimation of the accu- 
racy of the results given by the software, we have used 
SIMION8.0 to calculate the potential for infinitely long 
electrodes with hyperbolic sections in which case the po- 
tcntial can be exactly described by: <f> a (x, y) — x 9 f when 



Vq = and Uq — IV. Note that while the confinement is 
achieved dynamically, the anharmonic contributions to the 
potential do not depend on time and can be analysed as a 
static property of the trap. 

We can then calculate the average relative error Er(h), 
as a function of the node density h, using: 
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where <fi s (xi,yj) is the SIMION8.0 solution with Xi + i — 
Xi — y%+i — yi — h and n x , n y are the number of nodes. 

Figure [3] shows Er{h) for < x < r /2 and < y < 
ro/2 (the axes are oriented as in Fig[T]). Fig [3] also shows 
a fit using f(h) = a/h. The comparison of the two plots 
indicates that Er(h) closely follows the l/h law. The error 
bars correspond to one standard deviation, a. 

For the computation of the potential created by #1 
and #2, all possible symmetries have been used to reduce 
the required volume, as indicated in figure [2] by the dashed 
lines. This is important due to computer memory issues, 
which limit the mesh density which can be used for a given 
problem. In the present study, we were limited to a node 
density of 32 nodes/mm, for all computations carried out 
for a 3D volume. Figure [3] shows that this mesh density 
already provides an accuracy of 0.3 % for the hyperbolic 
trap. As the sizes of the studied traps are identical to 
the hyperbolic trap, the precision reached by SIMION8.0 
is of the same order. Therefore, the use of the potentials 
computed by SIMION8.0 in the study of the anharmonic 
content in a quadrupole linear trap is fully justified. 



1 In SIMION8.0, it is possible to define a 2D plane only, in which 
case, the program assumes this 2D section extends to [—00, +00] 
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Figure 3: (Color Online) Estimation of the error between 
the numerical solution given by SIMION8.0 and the known 
analytical potential in the case of infinite hyperbolic elec- 
trodes as a function of the mesh density h. The dots show 
E r (h) calculated by equation H whereas the curve plots 
f(h) — a/h, with a — 8.8 x 10~ 2 . The error bars corre- 
spond to lcr. 



2.2. Fitting the potentials 

An analytical two-dimensional solution of the Laplace 
equation for a quadrupole trap with circular sections ex- 
ists, if we assume infinitely long RF electrodes [TBI ]: 



(x,y) = Co + Real < C( 4m+2 )£ 



(4m+2) 



^ m— 



£ = x + iy 
(5) 

where C( 4m+2 ) are rea l coefficients and we have introduced 
a uniform contribution Co to equation ([5]) of fljj . 

As pointed out in the introduction, the term C§ van- 
ishes for infinitely long electrodes for r e = 1.14511ro flBT ]. 
To test our analysis of the potentials, we have computed 
the potential created by a 2D mass filter (infinitely long 
RF electrodes, without DC electrodes) for several values of 



the radius ratio r e /ro and found the C, 



(4m+2) 



coefficients 



of equation [5l by a fitting procedure analysed in the fol- 
lowing. The evolution of \Cq\ as a function of r e /ro shown 
in figure |4] shows that the Cq term can be decreased by 
3 orders of magnitude by choosing the appropriate radius 
ratio and confirms the zero crossing of Cq for a r e /ro value 
between 1.144 and 1.146. This accuracy corresponds to the 
maximal spatial resolution of our mesh for this problem (a 
high mesh density of 500 nodes/mm was possible due to 
the 2D nature of the problem). While this method does not 
provide a precise estimation of the "magic" radius ratio, it 
clearly agrees with the analytical result of r e = 1.14511ro, 
showing that SIMION8.0 produces a very good representa- 
tion of the potential created by an RF-trap. Higher order 
terms show a similar quantitative behavior as C§, with re- 
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Figure 4: (Color Online) Squares and right axis: absolute 
value of the first anharmonic coefficient Cg, versus r e /ro 
for ro = 4 mm. Blue circles and left axis: loss factor, L 
versus r e /r$. 



spective minima for Cio at r e /ro w 0.78, and r e /ro w 0.58 
for C14. 

In addition, figure S] shows the ratio between the ideal 
and the computed C2 coefficients which corresponds to the 
geometrical loss factor, C introduced in equation [TJ 
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For the chosen radius ratio, there is no loss of trapping 
efficiency induced by the circular shape of the electrodes 
(C ~ 1). We can observe that C takes values below 1 
for larger electrode radius, which for experimental rea- 
sons (decreasing inter-electrode distance) might not be 
an appropriate choice. In the following, all the studied 
trap configurations have the same radius tq — 4 mm and 
r e = 4.58 mm. 

The finite size of the rods and the presence of the DC 
electrodes introduce a dependency on z of the coefficients 
C(4 m+ 2). This dependency can be obtained by perform- 
ing a fit at each node along z to the numerical solution 
obtained by SIMION8.0. The trap parameters chosen are 
Vo =0, Uq = IV and Udc — 0. These values remain the 
same throughout the paper, unless specified otherwise. A 
truncated version of eqj5]is used to perform a least-squares 
fit [2(| • The points included in the fitting procedure belong 
to a 2D-square defined by < x < ro/2 and < y < ro/2, 
as r = ro/2 is considered as the maximum expected radius 
of the trapped ion cloud. 

To demonstrate the relevance of the fit, the x 2 of the 
fitting procedure for geometry #1.1 is plotted on figure [5] 
for each Zk and for increasing maximum order m max in- 
cluded in the fitting equation ([5])- These curves illustrate 
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Figure 5: (Color Online) Evolution of the \ 2 of the fit- 
ting procedure versus Zk/z , for increasing highest order 
m*max included in the fitting equation [5j The curves corre- 
sponding to m max — 2 and 3 overlap each other, showing 
that no gain is obtained by retaining orders higher than 



several features. First, they confirm that the Cq (m = 1) 
contribution does not modify the quality of the fit in the 
center of the trap, which is expected from the choice of 
the electrode radius r e that nulls the Cq contribution to 
the potential for infinitely long electrodes. Second, they 
demonstrate that adding higher orders than m = 2 does 
not lead to an improvement of the fit and in the following 
m = 2 is the highest order retained in the fitting func- 
tion. Furthermore, the x 2 parameter, which is lower than 
10 -9 for Zfe < zq/2 and m max — 2 increases drastically 
from Zk = zq/2 to zq. From Zk/zo = 0.6, one can see that 
adding higher order terms to the fitting equation does not 
result in an improvement of the \ 2 parameter. This means 
that equation ([5]) is not relevant in this part of the trap to 
reproduce the potential calculated by SIMION8.0 and that 
another analytical expression should be introduced to rep- 
resent the potential close to the DC electrodes. As we are 
interested in fitting the potential for z^/zq < 0.5, curves 
of figure [5] confirm the relevance of the fitting equation on 
the volume of interest. 

To compare the precision of the fitting procedure to 
the precision of the SIMION8.0 calculation, we define the 
relative error of the fit, Erfu(zk), taken as the relative 
difference averaged over the transverse plane, between the 
SIMION8.0 solution, (j) s (xi,yj, Zk), and the fitted poten- 
tial, 4> } (xi,y 0l z k ): 
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The relative error Erfu(zk) is calculated for the geome- 



try #1.1 and #2 for Zk < zq/2. For this relative error, 
the addition of the m — 2 contribution to the fitting po- 
tential allows reducing Erfn{zk) by a factor of 3 in the 
centre of the trap and by a factor of 2 at Zk = Zq/2, The 
relative error then reaches 5 x 10~ 6 in the center of the 
trap for both geometries and 8 x 10 -6 for geometry #2 
and 13 x 10 -6 for geometry #1.1 at Zk = Zq/2, This rel- 
ative error is three orders of magnitude smaller than the 
relative error one can expect from the SIMION8.0 calcula- 
tions. Nevertheless, fitting the calculated potential is rel- 
evant for comparison between different geometries as the 
relative precision of the fit is as good for the two geometries 
considered in the text and the precision of the SIMION8.0 
calculations should be identical as the basic geometry is 
the same for all the configurations compared here. Once 
all the coefficients Ci m +2{zk) arc determined at each Zk 
and for each configuration, it is possible to study the evo- 
lution of the anharmonic contribution along the trap axis 
and to compare the different configurations. 

3. Anharmonic contribution 

The relative anharmonic contribution, y, z) of the 
potential in the trap can be obtained using: 



with 
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where (j) r represents the real potential created by each ge- 
ometry and 4> q its harmonic contribution. A(f>f is the dif- 
ference between the SIMION8.0 solution and the fitted po- 
tential calculated using the C± m+ 2{zk) coefficients found 
by the fitting algorithm. A<p s represents the unknown dif- 
ference between the real and the SIMION8.0 solution, A£ 
is its relative value with respect to the harmonic contribu- 
tion (fig. 

In order to facilitate the visualisation of the behaviour 
of the anharmonic content, it is useful to calculate the 
averaged absolute value of £ along the z-axis: 



(C(**)> 
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The evolution of (((zk)) is shown in figure |6] for the 
configurations #1.1, #1.2 and #2. While all the configu- 
rations have comparable anharmonic content at the centre 
of the trap, strong differences appear between them when 
leaving the trap centre. As expected, due to the length of 
the RF electrodes with respect to the studied volume, con- 
figurations #1.1 and #2 show lower values of (C(zfc)) com- 
pared to #1.2. The increment observed from z^/zo > 0.5 
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Figure 6: (Color Online) Averaged anharmonic contribu- 
tion (C(zk)) along the trap axis for the configurations #1.1, 
#1.2, #2. Configurations #3 and #4 are introduced later 
in the text. Strong differences between the trap configura- 
tions appear from Zk/zo > 0.2. For clarity, only every 4 th 
node is plotted. 



Figure 7: (Color Online) Normalised potentials created 
by the DC electrodes on axis (0,0, z) for geometries #1 
(empty circles) and #2 (crosses) together with their re- 
spective fits to f{z) — az 2 (solid lines). Notice that the fit 
was performed considering only V(zk < 0.5zq). 



in the case of #1.1 can be explained by the unavoidable 
gap between RF and DC electrodes which generates a rise 
of the anharmonic content at the edge of the trap. The 
dip observed for (£#1.1(2*:)) is due to a zero crossing of the 
(^(xi,yj, Zk) in equation ((5]). In the area of interest (C( z fc)} 
remains rather constant along the trap axis for #2 and 
#1.1. 

Figure IH] shows that the truncation of the RF electrodes 
in configuration #1.2 leads to a substantial increase of 
the anharmonicity from Zk/zo > 0.3, which could limit 
the usable trapping volume. To demonstrate the effect of 
the boundary conditions, we have calculated the potential 
created by a configuration called #3, similar to #1.2 but 
with the DC-rods removed. We observe in fig [5] that 
this leads to a reduction of the the anharmonic terms, 
illustrating the importance of geometric aspects. 

We emphasize here that the only difference between 
#1.1 and #1.2 is in the voltage applied to the DC elec- 
trodes. Therefore, the discretization of the mesh by SIMION8. 
and so the calculation error is exactly the same for the 
two configurations. The fact that the expected behaviour 
of the anharmonic contribution is reproduced in each case 
re-enforces the validity of the assumption that SIMION8.0 
error contribution is the same for all the configurations 
and that the comparative study performed in this article 
is significant. 

4. DC-potentials 

In this section, we focus on the potential created by 
the DC electrodes inside the trap, as it also strongly de- 
pends on the chosen configuration. While at the very cen- 
tre of the trap, the potential can be approximated by an 



harmonic function, it is not necessarily true if larger vol- 
umes of the trap are considered. A harmonic potential 
shape can be desirable as it is often a necessary assumption 
when describing the ion dynamics or equilibrium proper- 
ties inside the trap, as for example to calculate the aspect 
ratio (radius over length) of an ion cloud [2l| . Figure [7] 



shows the potentials created only by the DC electrodes 
(Vb = Uq = 0) of the geometries #1 and #2, along the 
z-axis. For comparison, these potentials have been nor- 
malised using: 



VB° c rm (x,y,z) 



VDc(x,y,z)-V DC {0,0,0) 
V DC {0, 0,zq)- V dc (0,0,0) • 



(10) 



A fit to f(z) = az 2 using only the values of the potential 
in the interval < Zk/zQ < 0.5 (see fig[7j) clearly shows the 
divergence from an harmonic potential at the outer edges 
of the trap. In order to quantify the degree of harmonic- 
ity, we successively extended the fit, starting from the cen- 
Otre, until the correlation coefficient R 2 , commonly used to 
quantify the goodness of a fit, was worse than 0.9990. This 
was obtained at z^/zq — 0.240 and z^/zq — 0.597 for #1 
and #2 respectively. 

Another issue to take into consideration is the screen- 
ing of the DC-potential due to the presence of the RF elec- 
trodes, which is masked in figure [7] due to the normalisa- 
tion used. Assuming a potential depth AV^ = V(0, 0, z ) — 
V(0, 0, 0) of 1 V, we would need to apply to the DC elec- 
trodes 3.1 V for geometry #1, but 5.5 kV for #2. The 
requirement for very high voltages, without being criti- 
cal in many situations, can make option #2 inappropriate 
for an experimental setup where a high potential depth is 
needed or when designing an ion trap for space applica- 
tions where the weight of high power supplies eliminates 
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Figure 8: (Color Online) Description of the proposed al- 
ternative geometry, #4. The electric connections are the 
same as for trap #2, see table [T] for details. The DC elec- 
trodes geometry has been chosen so that the inner cut 
corresponds to the circumference inscribed in a square de- 
fined by the four rod centres. 

this configuration. On the other hand, #1.1 presents the 
disadvantage of the complexity of the electronics needed 
for driving the electrodes of the trap. 

The methodology presented so far, allows to study al- 
ternative designs and evaluate their performance. As a 
consequence, we propose an alternative geometry, which is 
introduced and analysed in the next section, showing that 
it is possible to achieve similar performances concerning 
low anharmonic contributions as #1.1 and #2, without 
the need for high voltages or advanced electronics. 

5. An alternative geometry 

An alternative geometry, denoted as #4 in the follow- 
ing, is described in fig [8] It uses continuous rods as in 
configuration #2 but with a different design of the DC 
electrodes that strongly decreases the screening of the RF 
electrodes. The evolution along the z-axis of the averaged 
anharmonic contribution (C(zfc)) is shown in Fig El We 
observe that the anharmonic content is comparable to the 
configuration #2 up to z^/zq ~ 0.4 and lower than #1.1 
for the rest of the trap. 

Regarding the DC-voltages, the correlation coefficient 
R 2 is found to be smaller than 0.9990 for z k /z Q < 0.28. 
Although the degree of harmonicity on the axial potential 
is not as high as for trap #2, this alternative geometry 
only needs an applied voltage of 37 V in order to have 
AV = 1 V. For some experiments, a smaller voltage can 
be preferable than a perfect harmonic behaviour on the 
axial direction. 

It is worth mentioning that the inner cut off of the DC 
electrodes corresponds to the circumference inscribed in a 
square centred on the electrodes as shown in figure [5] A 
geometry with the inner radius cut off equal to ro was also 
studied. The voltage required in that case for a potential 
depth Ay = 1 V was further reduced to 3.8 V. However, 
this reduced shielding is paid by an increase of the an- 
harmonicity which reaches the level of the worse geometry 



#1.2. This comparison shows that the anharmonic con- 
tributions are very sensitive to the position of the inner 
cut-off radius and configuration #4 is a good solution to 
reduce these contributions. 

6. Conclusion 

The anharmonic contribution to the potential created 
by a linear quadrupole trap has been studied for a large 
trap volume for several implementations of the DC and 
RF electrodes. We found that while the anharmonic con- 
tent is similar at the centre of the trap in all the different 
implementations, they behave extremely differently fur- 
ther away. The already operated configurations #1.1 and 
#2 as well as the proposed configuration #4 achieve sim- 
ilar good performances for Zk/zo < 0.5. Furthermore, for 
Zk/zo > 0.5 the RF potential of configuration #4 keeps a 
level of anharmonicity intermediate between #2 and #1.1. 
The rather complicated electronic set-up needed for the 
implementation of #1.1 and the need of high static volt- 
ages for the geometry #2, make these designs not suited 
for all experimental situations. The alternative geometry 
#4, proposed in this paper, represents therefore a trade-off 
between complexity of the implementation and low anhar- 
moncities, and should allow the envisaged trapping of a 
very large ion cloud. 
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